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Abstract 

This is the first of two papers where we discuss the limits imposed by 
competition to the biodiversity of species communities. In this first paper we 
study the coexistence of competing species at the fixed point of population 
dynamic equations. For many simple models, this imposes a limit on the 
width of the productivity distribution, which is more severe the more diverse 
the ecosystem is (Chesson, 1994). Here we review and generalize this analysis, 
beyond the "mean-field"-like approximation of the competition matrix used in 
previous works, and extend it to structured food webs. In all cases analysed, 
we obtain qualitatively similar relations between biodiversity and competition: 
the narrower the productivity distribution is, the more species can stably 
coexist. We discuss how this result, considered together with environmental 
fluctuations, limits the maximal biodiversity that a trophic level can host. 

1 Introduction 



One of the most striking characteristics of living systems is their amazing diversity. 
Theoretical ecologists have devoted much of their effort to explain why Nature is 
diverse and to identify the mechanisms that enhance or limit species coexistence. 
While the observations show that rich and diverse ecosystems are the rule and not 
the exception, stable, highly diverse systems rarely arise in mathematical models. 

In the early decades of theoretical ecology, emphasis has been placed in the 
principle of competitive exclusion. Stated qualitatively, this principle asserts that 
two or more species occupying the same ecological niche cannot stably coexist in the 
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same ecosystem. The first attempt to transform this principle into a mathematical 
theorem is due to Volterra (1928). Generalizations of Volterra's theorem to an S 
species ecosystem have been provided by several authors (MacArthur and Levins, 
1964; Rescigno and Richardson, 1965; Levin, 1970). It has been shown that S species 
cannot coexist at a fixed density if they are limited by less than S independent 
resources. 

This approach, however, presents two kinds of shortcomings. The first one is of 
mathematical nature. It has been shown that the theorem of competitive exclusion 
in the previous form does not hold if some conditions are relaxed. Most importantly, 
it no longer applies if the growth rates depend non-linearly on the resources and if 
coexistence at a fixed point is replaced by the more general condition of persistence 
(Armstrong and McGehee, 1980; Koch, 1974; Kaplan and Yorke, 1977). Such non- 
linear behaviour is to be expected in natural ecosystems. Similarly, the theorem of 
competitive exclusion is violated when age structured populations are considered, 
when resource is not uniform in quality (May, 1974; Diamond, 1978) or when physical 
space is considered (Sole et al, 1992). The prolonged coexistence of a large number 
of competitive plankton species has been justified through barriers to mixing of 
species in an otherwise homogeneous environment (Bracco et al, 2000). 

The second, more important difficulty is related to field observations, which can 
hardly ever claim to give a full account of all resources and the functional dependence 
of growth rates on them. Hence, the theorem of competitive exclusion has little 
actual predictive value. The simultaneous presence of similar species can always be 
ascribed to yet unknown resources or to unknown functional forms of growth rates. 

Field research has shown that coexistence of competing species is far from rare 
in real ecosystems. As a result, the conditions favoring coexistence are receiving 
increasing attention from the ecological community (McCann et al, 1998; Chesson, 
2000; Kokkoris et al, 2002; Roberts and Stone, 2004). These theoretical analysis 
and empirical observations suggest that the competitive exclusion principle could 
be replaced by a 'coexistence principle' (Boer, 1986). 

In the present suite of two papers, we analyze species coexistence for a class of 
models where competitive exclusion does not hold, since species growth is limited 
both by biotic resources (explicitly represented into the model) and by other limit- 
ing factors (implicitly considered), modeled as self-damping terms in the population 
dynamics equations. In the first paper we consider coexistence at the fixed point 
of population dynamics, and show how the combination of competition and envi- 
ronmental fluctuations limit the maximum amount of biodiversity that a trophic 
level can host. In the second paper we consider models of species assembly, kept 
away from the fixed point of population dynamics through the continuous arrival of 
new species due to immigration or speciation events. These models can reproduce 
species-area laws in good agreement with field observations. The coexistence condi- 
tion presented in this paper can be generalized to that situation as well. We define 
therefore an effective model of biodiversity across a food web, based on an approx- 
imation to the population dynamics equations and on the condition of maximum 
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biodiversity derived in this first paper. The effective model predicts that biodiversity 
as a function of the trophic level has a maximum at an intermediate level. 

The present paper is organized as follows. After introducing the population dy- 
namics equations that we use, we discuss in section 3 coexistence in a single trophic 
layer, in the framework of a mean-field approximation of the competition matrix. 
The coexistence condition that we derive imposing that all equilibrium densities 
are positive is equivalent to the one demonstrated by Chesson through the crite- 
rion of invasibility (Chesson, 1994; 2000). Several different models yield the same 
coexistence condition, be competition represented either explicitly in the popula- 
tion dynamics equations or implicitly, through the effect that other species have on 
resources. This condition states that species coexistence depends crucially on the 
distribution of rescaled net productivity, i.e. productivity not taking into account 
the competition load and rescaled through the carrying capacity: The productivity 
distribution has to get narrower for allowing the coexistence of a larger number of 
species. This result is also in agreement with theoretical studies stating that many 
species can coexist provided they are similar enough (Kokkoris et ai, 2002). 

The coexistence condition is then generalized to generic competition matrices, 
beyond the mean-field approximation. We show that the angle formed by the princi- 
pal eigenvector of the competition matrix, which we name the competition load, and 
the rescaled productivity must be narrow to permit species coexistence. This result 
allows to generalize the coexistence condition from one trophic layer to structured 
food webs, as discussed in Section 5. 

The coexistence condition alone does not impose a limit on biodiversity if the 
productivity distribution can be arbitrarily narrow. However, in natural ecosystems 
this distribution has a finite width due to unavoidable environmental fluctuations on 
time scales much shorter than those of population dynamics. This limits the maxi- 
mal biodiversity the system can host, and produces a typical shape of biodiversity 
versus trophic level (Lassig et ai, 2001), as we discuss in the companion paper. 

2 General framework 

We study here models of multispecies communities. A key ingredient in the models 
is biodiversity, meaning the number of reproductively separated populations in the 
environment. Biodiversity arises from a balance between the assembly process in- 
volving speciation and immigration and the extinction process driven by population 
dynamics, in a spirit similar to MacArthur and Wilson's theory of island biogeog- 
raphy (MacArthur and Wilson, 1967). Population dynamics is represented through 
generalized Lotka-Volterra equations of the form 
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where the community is formed by S interacting species, JVj represents the biomass 
density of species i, and N = {Ni, N 2 , ■ ■ ■ , N s }. 

The function ^(N) models prey-predator relationships, and is called Preda- 
tor Functional Response (PFR). It is zero for pairs of species not connected by a 
predator-prey relationship. If species % is a predator and species j is its prey, g^ is 
positive and represents the rate of prey consumption per unit of predator biomass. 
The sign is negative if species i is a prey and j is its predator. We assume in this 
case g ij (N)N i = —gji(N)Nj/rj. The term n < 1 is the efficiency of conversion of 
prey biomass into predator biomass. In the present paper we set r] = 1, in order to 
simplify formulas. In the companion paper the factor n will play a more important 
role, and it will be explicitly indicated. 

Several different mathematical forms of the PFR have been proposed and dis- 
cussed in the biological literature. We study here two cases: (i) prey dependent 
PFR linear in prey density, i.e. ^(N) = ^ijNj; (ii) ratio dependent PFR, i.e. g^ 
is a function of the ratio between prey and predator density (Arditi and Ginzburg, 
1989). 

The matrix models competition as a linear reduction in growth rates 1 . All 
its elements are non-negative, and the diagonal elements /3u are different from zero. 
These intraspecific competition terms play an essential role in allowing coexistence. 
Since competition for prey is already represented through the terms gij and the 
dynamics of prey species, the terms /3y stand for competition for resources not 
explicitly included in the model. Such terms naturally arise from "integrating out" 
some trophic links in a community. The term oii > accounts both for the death 
rate and for the energy consumption necessary for the activity of species i. 

Equations (1) are complemented with the condition that species below a critical 
density N c get extinct. This condition takes into account that species are composed 
of discrete entities and also mimics the effect of demographic stochasticity. 

We represent explicitly a single external resource, interpreted as abiotic and 
considered as an additional "species" N (t) (Caldarelli et al, 1998; Bastolla et al, 
2001). Its dynamics is chosen in such a way that species feeding on it feel an indirect 
competition effect. The qualitative behavior of the model does not depend on the 
detailed dynamics of the abiotic resource, as far as competition in the first trophic 
level is properly represented. The external resource introduces a new scale of density 
R. The dimensionless quantity R/N c ^$> 1 is an important control parameter in the 
system. 

lr This linear model can be understood as the linearized version, around an equilibrium point, of 
a more general model, as far as equilibrium properties are captured in the linearized representation. 
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3 One-layer competition and productivity distri- 
bution 



In this section, we reformulate results showing that the coexistence of several com- 
peting species tends to equalize their net productivity: the more coexisting species, 
the more similar their productivity should be (Chesson, 1994; 2000; Lassig et al., 
2001; Kokkoris et al, 2002). This is done here imposing that the S competing 
species coexist at a fixed point of population dynamics, with all densities positive 
and larger than the threshold for extinction, N c . We adopt a simple mean-field 
approximation of the competition matrix, which will be relaxed in next section. 

The condition derived in this way is qualitatively equivalent to a condition de- 
rived through the more general requirement of invasibility (Armstrong and McGehee, 
1980; Chesson, 1994; 2000). Moreover, modelling competition through explicit in- 
teraction terms or implicitly, through the dynamics of the common resources, leads 
to the same condition. 

3.1 Direct competition: mean-field approximation 

Our study of competition starts from the simplest model 



The quantity Pi is assumed to be independent of species density, and represents the 
intrinsic growth rate of species % in the absence of competition. 

We assume that the competition matrix is symmetric and all its elements 
are non-negative. We further assume for convenience that the matrix is positive 
definite. It can be proven that these hypothesis hold if the competition terms arise 
in an effective way through the dynamics of underlying resources. The stability 
properties of the fixed point N* = Y,k(P)ik Pk are governed by the community 
matrix (May, 1974), which in the present case is = —N*^. For positive definite 
competition matrices, positivity of all the iVj implies that the community matrix 
is negative definite, thus the system is locally stable. Furthermore, under these 
conditions, it is possible to construct a Ljapunov function (MacArthur and Levins, 
1964; May, 1974), which guarantees that the equilibrium point is globally stable as 
well. Thus, we can ignore coexistence along periodic as well as chaotic orbits. 

We parameterize the competition matrix as 



where pij G [0, 1] is a dimensionless quantity that we call ecological overlap (or niche 
overlap) and describes the similarity in the use of resources between species i and 
j. One has clearly pa = 1. 



1 dNj 
Ni dt 
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The parameters flu can be absorbed introducing the rescaled variables rii = 
{PuY^Ni, pi = {f3a)~ l l 2 Pi. The variables pf have dimensions of biomass per area 
per time, the same dimensions of productivity. By analogy and as a shortening, we 
call the pj's productivities instead of rescaled growt rates, but they should not be 
mistaken for productivities measured in field studies. In terms of the new variables, 
the fixed point equations have the form 

Pi = J2 Pn n i ■ ( 4 ) 

3 

We start by considering a mean field approximation of the competition matrix, 
addressing the general case in the next section. The mean field approximation 
consists in assuming that all non-diagonal elements are equal: pa = 1, p^- = p for 
i 7^ j, 2 such that the solution of the fixed point equations reads 

1 ( (p) \ ^ 

ni "(i-p )r i + a - poVpos) ' [b) 

where ( ) indicates the average over the S species in the community. All equilibrium 
densities are positive and above the threshold n c if and only if 

(P) ~ Pi < 1 - nj (n) 

(p) - 1 + Sp /(1 - p ) ' {) 

For n c = 0, this result is equivalent to the condition derived by Chesson (1994) 
imposing invasibility of the system. 

The above condition only affects explicitly the pj's smaller than the average, 
but it is easy to see that it implies a condition on the variance of the productivity 
distribution. In fact, multiplying both sides by the quantity (1 — (pi — (p))/(p)), 
which is always positive if the variance of the distribution is small enough, and 
averaging over all species, we find 

(P 2 ) - (P) 2 < (1-Po) (7) 



(p) 2 -Spo + 0—po) 

In other words, the coexistence condition requires that all p^s are very close to 
the average value when the number of species S is large. The maximal negative 
difference from the average is of order 1/S and the standard deviation is of order 
l/y/S. Notice that, if the mean overlap p equals unity, all p, L must be identical. 
This result corresponds to the formulation of the theorem of competitive exclusion 
in this framework. 

For a resource rich system where the average reduced density (n) is well above 
the extinction threshold, the quantity So = (1 — Po)/po defines an intrinsic scale of 
biodiversity at which the productivity distribution gets pretty narrow. For much 

2 In the context of May and MacArthur's theory of competition, this hypothesis is realized when 
the niche space has many dimensions (May and MacArthur, 1972). 
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larger number of species S, the average reduced density (n) = (p) / (1 + (S — l)po) ^ 
(p)/(l + S/So) is close to the extinction threshold, thus the r.h.s. of Eq. (6) becomes 
very small and the condition of coexistence becomes rather stringent. Therefore, the 
effect of positive n c on biodiversity is only important for very diverse ecosystems. 

3.2 Resource mediated competition, prey dependent PFR 

We now consider S basal species feeding on a single external resource N with a 
linear, prey dependent PFR: gij({Nk}) = IjNj. Competition is induced through the 
dynamics of the common resource N and through the limiting factors not explicitly 
considered in the model. The population dynamics equations are 

3 

S 

RPo-PoN -J2liNi. (8) 
i=i 

Different equations for Nq give qualitatively similar results, as far as iVo is consumed 
by all the competing species. We consider the fixed point equations and substitute 
for N , finding 

^-Z(^ + ^)^ = 0, (9) 

where we have defined the rescaled variables 7- = •yi/y/Pu, a[ = ai/y/Pu, pi = 
7j'i? — a[ and rii = iVj-^/Ai- The fixed point equations are thus equivalent to the 
equations obtained by considering direct competition, with an effective competition 
matrix that is symmetric and positive definite if j3ij is such. 

We now adopt the mean- field approximation of the previous section, pa = 1, 
Pij — Po < 1 (i 7^ j)- The approximation is applied only to the matrix p^-, since 
the other part of the competition matrix depends on the productivity vector, nev- 
ertheless the result turns out to be equivalent to the previous one. Substituting 
li = (Pi + ol')/R in Eq. (9), we get 



1 dNj 
Nt dt 



iVn dt 



(1 - p )nj = pi(l - S(n)a) - aiS(n)a - p S(n) , (10) 

where a = (y'n)/ (n)(3 R. We then neglect the dependence of a on (n) and solve 
for this variable. Rearranging the various terms, we find again that the coexistence 
of S species is only possible if the minimal productivity differs from the average at 
most in an amount (1 — p)j ' pS: 

Pi - (P) fljgO <-(<*') l-n c /(n) 

(p) a{a') + p + (l-p )/S (a') ~ 1 + Sp/(1 - p) ' 1 ' 
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where the effective competition overlap p is defined as 



f±^M<i (12) 

1 + a {a') 

This coincides with Eq. (6) in the limit a(a') = (a') ('j'n) / (n) (3qR — > 0, in which 
resources are not limiting. In the opposite limit the condition on minus the energy 
consumption rate — a[ becomes more demanding than the condition on the p^'s, and 
p tends to one. Notice that the competition overlap p is different from zero even if 
the explicit competition overlap po vanishes. 

3.3 Resource mediated competition, ratio dependent PFR 

We now assume that the PFR depends on the ratio between the density of the prey 
and the density of its predators (Arditi and Ginzburg, 1989) as 

b CiN 

(see also Schreiber and Gutierrez, 1998, Sole et al, 2000, Drossel et al, 2001). In 
this expression species "0" is the prey and the sum in the denominator runs over its 
predators, represented as the set -P(O). In the case of just one predator, Eq. (13) can 
be seen as a Holling type III PFR (Holling, 1959), where the prey density at which 
the functional response saturates is proportional to the predator's density. Notice 
that, unlike prey dependent PFR, the ratio dependent PFR does not contain any 
externally specified biomass scale. For this reason, it is in some sense simpler than 
prey dependent PFR (Bastolla et al, 2001). We do not aim at discussing which form 
of the PFR is most suitable to explain observational data (Abrams and Ginzburg, 
2000), but just use it as a second example of population dynamics where, aside of 
ratio dependence, competition for preys is explicitly represented through the sum 
over predators in the denominator. Since competition is now explicitly represented, 
we don't need an additional equation for the external resource N and we assume 
that it renews rapidly enough such that its density remains constant. The dynamical 
equations have now the form 

1 dNi b Cl N _ 



N t dt bNo + Y.jCjNj 

N (t) = R. (14) 

To simplify formulas, we use the rescaled variables c\ = Cij '^/Pu, rii = iVjA/Ai and 
apply the mean field approximation to the matrix fy. The solution of the fixed 
point equations is then 

1 c' - (c') 

n* = (n) + - . n Ct //. , (15) 
1 w (1 - po) 1 + Sa{n) v 1 
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where now a = (c'n) / (n)bR. To simplify formulas, we consider the case of a pro- 
ductive system where the variation in the a^s can be neglected. The condition that 
all species are viable yields 




< 



l + Sp/(l-p) ' 



1 - nj (n) 



(16) 



where 



P = 



p + (a')/PS(n) 
1 + (a')/pS(n) 



(17) 



is the effective competition coefficient and S(n) can be obtained solving a second 
order equation, whose root remains finite in the limit S — > oo. 

Equations (6), (11), and (16) are equivalent coexistence conditions for three 
different population dynamics models. They all show the quantitative dependence 
between the distribution of productivities and biodiversity. Qualitatively identical 
coexistence conditions were obtained through the invasibility criteria, i.e. imposing 
that an invading species (with very low density) has a positive growth rate (Chesson, 



4 One layer competition without mean-field 

In this section, we analyze the coexistence condition without relying on the mean 
field approximation of the competition matrix. Our starting point are the fixed 
point equations (4), where rescaled variables are used. To simplify the presentation, 
most calculations are reported in Appendix A. Here we only show and discuss the 
final results. 

We will use the spectral representation of the matrix p^, which we assume to 
be symmetric and positive definite. Therefore, all of its eigenvalues \ a are real 
and positive, and its S eigenvectors u a form an orthonormal system. In graph 
theory (Bollobas, 1998) the matrix p^ is called the adjacency matrix. Its maximal 
eigenvalue is defined by the property Ai J2i v f > Pij^iVj for every vector v. The 
equality holds if and only if v is proportional to the principal eigenvector u 1 . It holds 
that Ai > Pij/ S, so that the main eigenvalue is expected to be proportional to the 
number of species S. Moreover, all components of the principal eigenvector either 
have the same sign (we can choose it arbitrarily to be the positive sign) or are zero. 
We will assume that the graph cannot be separated in disconnected components. 
If this is not the case, the analysis can be applied separately to each disconnected 
component. This hypothesis implies that all components of the principal eigenvector 
are strictly positive: u\ > 0. We call the principal eigenvector of the competition 
matrix the competition load. 

It is useful to define the average of the 5 — 1 eigenvalues excluding the principal 



2000). 
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one, which we will denote with the symbol 1 — po : 



a=2 



It follows from the condition on the trace of the matrix pij that J2 a = S, so 
the principal eigenvalue Ai can be expressed as 

Ax = S Po + (1 - p ) • (19) 

We now use for simplicity the notation v a = J2j VjU^ for the projection of a vector 
v along the a-th eigenvector, u a . We show in Appendix A that the coexistence 
condition (6) can be written, for a general competition matrix, as 

gj£LlW! < A. (l-^) (20) 

(P 1 ) 2 - Sp„ + (1- P(1 ) 1, (n)j • y ' 

where A2 is the second largest eigenvalue, which is of order one. 

If the matrix p^- has the mean field form, Eq. (20) is identical to Eq. (7). In 
fact, in the mean field case, it holds that X a — (1 — po) for a > 2. Moreover, in 
this symmetric system the principal eigenvector is uniform, u\ = 1/V~S, so that 
p x = y/S(p), yielding Eq. (7). 

Eq. (20) tells us that the productivity vector, p iy has to be almost parallel to 
the competition vector u\ in order to allow coexistence of all species: Species with 
larger competition component uj need a comparatively larger productivity to sur- 
vive. The angle between the productivity vector and the competition vector has to 
become narrower as competition becomes stronger, until only a perfect coincidence 
guarantees the survival of all species. The strength of competition increases with 
the overlap parameter 1 — A2, with the ratio n c / (n), and with the largest eigenvalue 
Ai ~ Spo- Therefore, competition becomes more severe as the number of species 
increases, unless the competition matrix is sparse. Since the resources on which 
the S species feed are finite, we expect on ecological grounds that the strength of 
the competition increases as more species are packed in the ecosystem. This can be 
shown in some mechanistic models of competition, in which the matrix is derived 
from the underlying dynamics of the resources. 



5 Competitive coexistence in a food web 

We now turn to the more general case of an L— levels food web. After integrating 
out the upper and lower levels, the effective equation for the fixed point density of 
species at a given level has the familiar form of a competition equation, and displays 
the same qualitative behaviour discussed for the case of the single level. 

As in our previous work (Lassig et ai, 2001), we assume a hierarchical trophic 
organization, so that species at level / feed only on species at level I — 1 and compete 
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only with species at their same level. Species at level L are top predators. Level 
zero can be interpreted either as basal species or as abiotic resources described by 
an effective equation. The dynamical equations, using linear, prey dependent PFR, 
read 



N (l) dt 



ETg^-a? (21) 



where the superindex stands for the level where the species belongs. The matrix 7® 
represents predation from level / to level I — I. It vanishes identically for I — and 
I — L + 1. The parameters af 1 stand for energy consumption and death rates of 
species at level / for I > 1, whereas — is the fictitious growth rate of the resource 
at level zero. The direct competition for species at level / is represented through the 
matrix ft®. As above, we assume that it is symmetric and positive definite. 

We now consider the fixed point equations. For every level /, we solve for the 
densitites of species at different levels I' ^ I and substitute, thus getting fixed point 
equations for species at level I with the form of effective competition equations with 
a symmetric competition matrix: 

E4°< } = p P- (22) 

j 

Three different integration schemes yield different equations of the same form (22). 
In the scheme described in Appendix B, densities at levels V < I are recursively 
solved starting from level zero, and densities at levels I' > I are recursively solved 
starting from the maximum level L. Levels above and below I (predators and preys) 
contribute both to the competition matrix and to the productivity vector. Therefore, 
as it is known, competition can be induced not only through common preys, but 
also through common predators. 

In a second possible scheme, levels above I only contibute to as an effective 
energy consumption term, whereas levels below / contribute both to the competition 
matrix and to the productivity vector as a growth term. Finally, in the third scheme, 
levels above / contribute both to the effective competition and to the productivity 
vector, whereas species at levels below I only contribute to the productivity. In all 
three cases, the effective equations for species at level / have the form (2) of effective 
competition equations, where the competition matrix and the productivity vector 
depend on the properties of species at the other levels. 

The result of the previous section implies that also for a structured food web the 
effective productivity vector at level / must form a narrow angle with the principal 
eigenvector of the effective competition matrix. 
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6 Environmental variability and biodiversity 



The coexistence conditions derived above yield a natural scale for assessing the 
degree of biodiversity of an ecosystem 3 . Ecosystems with S <C (1 — p)/p are loosely 
packed. They impose mild conditions on the productivity distribution, and can 
easily incorporate new species. Ecosystems with S ^> (1 — p)/p are tightly packed. 
There, the productivity distribution is subject to strict constrains and incorporation 
of new species is very difficult. Which ecological mechanisms distinguish between 
the two kind of situations? On the one hand, as we will discuss in a next paper, the 
immigration and speciation rates play a key role in smoothing the effects of condition 
(6) and allowing a large number of species into the system. On the other hand, an 
important role in determining biodiversity is played by environmental fluctuations 
(Hutchinson, 1961; Chesson, 2003a). 

The growth rates considered in the previous models are subject in real systems 
to environmental fluctuations, such as for instance variation in rain levels, temper- 
ature and daily light, or forest fires for plant communities, natural obstacles to the 
movement of animal species, or fluctuations in the number of interacting popula- 
tions also affected by the environmental noise. Some of these variables fluctuate over 
time scales much shorter than the characteristic time scales of population dynamics. 
We thus interpret the growth rates p^s as time averaged quantities with superim- 
posed fast fluctuations due to environmental variability. These preclude that the 
productivities take identical values, even for populations with identical ecological 
properties. 

To take into account these fast environmental fluctuations, we assume that they 
tend to broaden the productivity distribution so that there is a minimum width 
given by the equation 



The factor 1 — 1/ S ensures that the condition (23) is satisfied for S — 1, when P m i n 
and (P) coincide. The term < A < 1 is an effective measure of environmental 
variability. For A = the condition (23) is always satisfied, while for A = 1 it 
imposes a very small P min (recall the P's are positive quantities). Values of A > 1 
can not be realized in the large S limit. 

Combining Eq. (23) with the coexistence condition (6), we get a limit to the 
maximal biodiversity that the ecosystem can host. In fact, Eq. (6) requires that 
the productivity distribution gets narrower as biodiversity increases, but at some 
point the largest difference tolerated becomes of the order of the minimal difference 
compatible with the actual environmental variability. The maximal biodiversity is 

3 This scale is relevant in ecosystems where (N) ^> N c , otherwise biodiversity is controlled by the 
threshold density. We are assuming here that our ecosystem is large enough so that the condition 
(N) 3> N c holds where the number of species is of the order (1 — p)j p at which competition effects 
are important. 




(23) 
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a function of the competition overlap p, of the environmental variability A, and of 
the ratio between the average biomass and the threshold for extinction, n c j (n) , and 
is given by the inequality 




(24) 



leading to a second order inequality that can be simplified for large S as 



S < 1 + 



l-p \ ( l-A-n c /{n) 



T) V a 



) 



(25) 



For large competitive overlap p close to unity or large variability A ~ 1, only one 
species can survive in the long run. If both the overlap and the variability are small, 
on the other hand, the maximal number of species can be rather large. Only in case 
of small variability the ecosystem can be tightly packed. 

7 Discussion 

In this work, we have considered simple models of competition, either represented 
as explicit terms in the population dynamics equations or effectively introduced 
through the dynamics of common prey and predator species. We have focused 
our attention on the fixed points of population dynamics, imposing the condition 
that all species at equilibrium have positive densities. For the simplest model of 
competition, this condition ensures that the fixed point is both locally and globally 
stable, but for more complex situations this represents an oversimplification. We 
will consider in the companion paper systems far from equilibrium, and argue that a 
suitable modification of the coexistence condition derived at the fixed point remains 
generally valid. 

We first analyzed coexistence at the fixed point, adopting a mean field approx- 
imation of the competition matrix. This approximation is expected to hold when 
there is a main nutrient on which all species are dependent, or when the niche space 
has many dimensions. 

Within the mean-field approximation, coexistence of S species competing with 
each other implies a condition on the distribution of their rescaled growth rates, 
Pi = Pi/ yffiii, where flu is the inverse of the carrying capacity: The width of the 
distribution has to shrink with increasing number of species S, competition overlap 
p or ratio N c / (N) of the minimal viable density to the average density. Notice that, 
if the carrying capacities vary from one species to another, for instance depending 
on body size, this conclusion does not apply to growth rates prior to rescaling. 

This coexistence condition is equivalent to the one obtained by Chesson (1994; 
2000) by imposing the requirement of invasibility, which demonstrates its robust- 
ness. We found formally identical coexistence conditions in models with explicit 
competition terms and in models where competition is induced by the dynamics of 
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shared resources, within a single trophic layer and within a whole structured food 
web. 

The mean field approximation is not appropriate to describe more complex com- 
petition matrices, for instance when there is a continuum of resources distributed 
along one dimension (May and MacArthur, 1972), or when species can be grouped 
according to their ecological similarity. We have therefore generalized the coexis- 
tence condition beyond mean-field competition matrices. In the general case, an 
important role is played by the principal eigenvector of the competition matrix, 
which we named the competition load. The angle between the vector of competi- 
tion loads and the vector of rescaled growth rates has to be narrower for stronger 
competition, i.e. for increasing average competition overlap and number of species. 
This generalizes the condition on the variance of the distribution of rescaled growth 
rates to the case of a generic matrix. 

The coexistence condition Eq. (7), stating that the variance of the rescaled 
growth rates must decrease as 1/S due to competition, is reminiscent of May's 
theorem, according to which the coexistence of S species randomly interacting re- 
quires that the variance of their interactions vanishes as 1/S (May, 1974). In fact, 
May's theorem was derived imposing local stability of the fixed point of population 
dynamics and, for the simplest models of competition, the condition that all densi- 
ties are positive at the fixed point also implies local stability. The condition for the 
coexistence of species whose rescaled growth rate lies below the average is stricter: 
they must differ from the average at most by an amount of order 1/S, which is 
smaller than that imposed by the condition on the variance. This appears more 
demanding than the condition in May's theorem, probably because here all species 
are in competition with each other instead of being randomly interacting. 

The coexistence condition that we described does not limit biodiversity, as long 
as the angle between rescaled growth rates and competition loads can be made 
arbitrarily small. However, environmental fluctuations, due to biological or abiotic 
processes with time scales shorter than those of population dynamics, necessarily 
limit the possibility to fine tune ecological parameters in order to accomodate a 
larger biodiversity. We include environmental variations in our competition model 
in an effective way, as a force broadening the distribution of rescaled growth rates. 
Considering this new ingredient sets a limit on the maximal biodiversity that a 
system of competing species can host. In the companion paper, we will combine 
this limitation to "horizontal" biodiversity imposed by competition, with either 
dissipation of energy or growth of perturbations in the "vertical "direction along the 
food chain. 

Environmental changes characterized by longer time scales, as seasonal changes 
or the storage effect, have not been included in our analysis, although they certainly 
affect coexistence conditions. Interestingly, in some cases the adaptative responses 
of species to these environmental fluctuations with time scale comparable to that of 
population dynamics are predicted to enhance species coexistence (Chesson, 2003a; 
2003b). This constrasts with the prediction presented here that a fast fluctuating 
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environment would set a limit on biodiversity. It would be therefore quite desir- 
able to set up a general theory of how environmental noise on different time scales 
modulates biodiversity. 
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Appendix A: Coexistence condition for a generic 
matrix 

Our starting point is Eq. (4), = £) ■ PijUj, and we want to demonstrate the coex- 
istence condition (20) using the spectral properties of the competition matrix p^. 
For simplicity of notation we shall indicate with a superscript a the component of a 
vector in the direction of the a-th eigenvector of the matrix p^-, n a = J2j n j u "j- A a 
denotes the corresponding eigenvalue. Using this notation, Eq. (4) can be written 
in the eigenvector basis as p a = \ a n a , from where we get the equilibrium biomasses 

n, = £^>n c . (26) 

a=l A " 

We now bring to the r.h.s. the contribution of the principal eigenvector u\ and multi- 
ply both sides times the quantity \{j>i — p 1 u})/p 1 — uj], which, as we will demonstrate 
at the end, is always negative if system size is large enough. We get 

^ P a uf f pi- p^u] \ ( p l u\\ ( pi-p x u\ \ 

k A. I fuj V " 1^ Ai ) [ p^uj V • {27) 

Summing over i, exploting the ort honor mality of the eigenvectors, and rearrange 
the factors, we obtain 

E ^ < ^ - n c (2p V) - <p>) • (28) 

Here brackets denote average over the S species. It holds (p) = Y, a P a ( ua )- Ne- 
glecting , just for the sake of simplifying the final formula, the contribution of the 
eigenvectors with a > 2 to the mean productivity, which is small since the compo- 
nents p a and the mean eigenvectors (u a ) are much smaller than the corresponding 
quantities for the principal eigenvector, we get 
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Noticing that A 2 is the largest of all eigenvalues with a > 2, and that I]f =2 (p a ) 2 = 
Y,i=iPi — (p 1 ) 2 , and subsituting the expression for Ai, we finally get Eq.(20): 

This expression demonstrates our previous statement that the quantity (pi—p 1 u})/p 1 u 1 
is always smaller than one if system size is large enough, which is the property that 
we have used for deducing the above formulas. 



Appendix B: Effective competition in a food web 

Starting from the food web equations (21), we want to obtain effective fixed point 
equations for level / in the form (22). First, we write the fixed point equations in 
matrix notation: 

( 7 (D) N d-i) _ a (D _ (0(D) N (0 _ ( 7 «+D) T N ( ;+1 ) = . (31) 

Bold face symbols represent column vectors, while the other symbols represent 
matrices. An upper T indicates tranposition of a matrix. The boundary conditions 
are 7 ( L+1 ) = (the top level L has no predators) and a formally identical equation 
with constant growth rate for the renewable abiotic resources at level 0: 

7 (o) jR _( 7 (D) T N (i) = 0(0)^(0) _ (32) 

We now focus on an intermediate level /. Species below level I can be integrated 
out solving for their equilibrium densities starting from level upwards. Species 
above level / can be integrated out similarly starting from the maximum level L. 
Both kinds of species contribute to the effective competition matrix and effec- 
tive productivity vector of Eq.(22). The recursive equations to obtain these 
quantities are 



C« = ( 7 «)(m('- 1 ))" 1 ( 7 ») T + / 3(0 

+ ( 7 (m) ) T (M (m) ) _1 ( 7 (m) ) (33) 

p(0 = (^(M^'V^-c^ 

+ (/+ 1 )) T (M( i+1 )) _1 Q^ 1 ) (34) 



Quantities with the tilde are recursively obtained from the upper levels, starting 
from level L: 
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m = + ( 7 (m) ) T (M^y 1 ( 7 { ' +i) ) (35) 

Q« = -«<') + ( 7 { ' +1) ) T (M (m) ) _1 Q(' +1 ) . (36) 

The boundary condition is (j&( L ^ = (3^ L \ = —a^ L \ Quantities without the 
tilde are similarly obtained starting the recursion from level zero: 

M« = ( 7 «) (M^y 1 ( 7 «) T + /?« (37) 
Qf) = ( 7 f))(Mf- 1 ))" 1 Qf- 1 J-a« (38) 

The boundary conditions are in this case (^M^^j = (3^ and Q(°) = 7 ^i?. 

The effective competition matrix is positive defined if all of the (3 and 7 
matrices are so. 
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